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A KINETIC  APPROACH  TO  COLLISION  PROCESSES  IN  GASES. 


I.  SMALL  AMPLITUDE  PROCESSES  IN  CHARGED  AND  NEUTRAL 

ONE  COMPONENT  SYSTEMS 

by 

* **  * 

P.  L.  Bhatnagar,  E.  P.  Gross  and  M.  Krook 

Abstract:  A kinetic  theory  approach  to  collision  processes  in  ionized  and  neutral 
gases  is  presented.  This  approach  is  adequate  for  the  unified  treatment  of 
the  dynamic  properties  of  gases  over  a continuous  range  of  pressures  from 
the  Knudsen  limit  to  the  high-pressure  limit  where  the  aerodynamic  equa- 
tions are  valid.  It  is  also  possible  to  satisfy  the  correct  microscopic 
boundary  conditions.  The  method  consists  in  altering  the  collision  terms  in 
the  Boltzmann  equation.  The  modified  collision  terms  are  constructed  so 
that  each  collision  conserves  particle  number,  momentum  and  energy;  other 
characteristics  such  as  persistence  of  velocities  and  angular  dependence 
may  be  included.  The  present  article  illustrates  the  technique  for  a simple 
model  involving  the  assumption  of  a collision  time  independent  of  velocity; 
this  model  is  applied  to  the  study  of  small  amplitude  oscillations  of  one- 
component  ionized  and  neutral  gases.  The  initial  value  problem  for  un- 
bounded space  is  solved  by  performing  a Fourier  transformation  on  the  space 
variables  and  a Laplace  transformation  on  the  time  variable.  For  uncharged 
gases  there  results  the  correct  adiabatic  limiting  law  for  sound-wave  propa- 
gation at  high  pressures  and,  in  addition,  one  obtains  a theory  of  absorption 
and  dispersion  of  sound  for  arbitrary  pressures.  For  ionized  gases 
* Astronomy  Department,  Harvard  University,  Cambridge,  Massachusetts. 
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the  difference  in  the  nature  of  the  organization  in  the  low-pressure  plasma 
oscillations  and  in  high-pressure  sound-type  oscillations  is  studied.  Two 
important  cases  are  distinguished.  If  the  wavelength  of  the  oscillations  is 
long  compared  to  either  the  Debye  length  or  the  mean  free  path,  a small 
change  in  frequency  is  obtained  as  the  collision  frequency  varies  from 
zero  to  infinity.  The  accompanying  absorption  is  small;  it  reaches  its 
maximum  value  when  the  collision  frequency  equals  the  plasma  frequency. 
The  second  case  refers  to  waves  shorter  than  both  the  Debye  length  and 
the  mean  free  path;  these  waves  are  characterized  by  a very  heavy  ab- 
sorption. 

INTRODUCTION 

The  dynamical  theory  of  gases  may  be  studied  from  two  points  of  view. 
One  may  take  as  starting  point  the  macroscopic  equations  of  aerodynamics  with 
the  density  p , mass  velocity  q and  temperature  T as  independent  variables, 
and  involving  various  coefficients,  e.g.,  viscosity,  heat  conduction,  etc.  On 
the  other  hand,  one  may  use  a more  fundamental  and  general  microscopic 
formalism.  The  most  fruitful  of  such  formalisms  available  at  present  is  that 
in  terms  of  one  particle  distribution  functions  satisfying  integro-differential 
equations  of  the  Boltzmann  type. 

The  aerodynamic  equations  are  adequate  for  treating  a certain  wide 
range  of  problems  in  gas  dynamics  and  most  of  the  developments  of  the  theory 
have  been  in  this  domain.  In  this  range  the  microscopic  theory  would  not 
yield  significantly  different  results.  In  fact,  subject  to  certain  conditions  which 
delimit  their  range  of  applicability,  the  aerodynamic  equations,  together  with 
explicit  formulae  for  the  various  coefficients  entering  into  them,  are  derivable 
from  the  microscopic  theory.  A detailed  account  of  this  derivation  is  presented 
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in  the  book  of  Chapman  and  Cowling.^ 

There  are,  however  many  important  situations  in  which  the  macroscopic 
theory  does  not  give  a correct  description.  During  recent  years  interest  in 
such  problems  has  increased  considerably.  An  extreme  example  of  the  breakdown 
of  the  aerodynamic  theory,  viz. , in  the  Knudsen  region,  has  been  known  for  some 
time;  here  the  mean  free  path  of  a molecule  is  large  compared  to  some  linear 
dimension  of  the  apparatus  and  behavior  at  the  boundaries  becomes  important. 
Other  examples  are  provided  by  high-frequency  sound  waves  in  a rarefied  gas 
and  plasma  oscillations. 

The  solution  of  the  Boltzmann  equation  is,  in  general,  a matter  of  consider- 
able difficulty  even  in  cases  corresponding  to  the  physically  simplest  situations. 
Significant  progress  has  been  confined  practically  to  the  study  of  two  limiting 
cases  in  which  two  different  approximation  procedures  can  be  applied.  A cri- 
terion for  the  range  of  validity  of  the  approximate  methods  is  provided  by  the 
comparison  of  some  characteristic  time  r(or  characteristic  length  L.)  with  the 
average  time  TQ  (or  mean  free  path  L.c)  between  molecular  collisions. 

For  high  density  ( T » TQ  or  L >>  L<c)the  Enskog-Chapman  theory  may 
be  used.  The  first  approximation  of  the  theory  consists  in  assuming  local  thermo- 
dynamic equilibrium  and  a common  drift  velocity  for  all  molecular  species.  The 
next  approximation  corrects  the  distribution  function  by  terms  proportional  to 
the  first  derivatives  of  temperature,  velocity  and  density;  this  corresponds  to 
the  aerodynamic  equations  with  coefficients  of  heat  conduction,  viscosity  and 
diffusion.  The  high  density  region  (r  » T^.)  is  in  fact  the  range  in  which  the 
aerodynamic  equations  provide  an  adequate  description.  Higher  approximations 
of  the  E-C  theory  lead  to  correction  terms  proportional  to  higher  derivatives  of 

1)  S.  Chapman  and  T.  G.  Cowling,  "The  Mathematical  Theory  of  Non-Uniform 
Gases",  Cambridge  University  Press,  London,  2nd  Edition,  1952. 
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T,  q,  p . The  successive  approximations  of  the  E-C  theory  correspond  to  an 
expansion  of  the  distribution  function  in  powers  of  the  mean  free  path  L£.  If 
we  consider  sound  waves  with  wavelength  !_.=>>  Lc>  the  first  and  second  approxi- 
mations are  already  sufficient  to  give  all  significant  features  of  the  process. 
However,  when  L becomes  comparable  with  Lc>  it  is  necessary  to  go  to  the 
third  and  even  higher  approximations  to  obtain  adequate  results;  as  discussed 
by  Kerzfeld  ^ important  dynamic  contributions  appear  only  in  approximations 
higher  than  the  second.  The  third  approximation  already  involves  formidable 
labor  and  has  been  used  to  solve  only  the  simplest  problems.  Consideration  of 
higher  approximations  is,  in  any  case,  of  doubtful  value  as  the  entire  procedure 
breaks  down  in  just  the  range  where  the  contributions  from  these  higher-order 
terms  become  important.  In  addition,  the  boundary  conditions  in  many  problems 
cannot  be  specified  properly  within  the  scope  of  the  E-C  formalism.  A different 
approach,  using  expansions  in  terms  of  Hermite  polynomials  in  velocity  space, 
has  been  given  by  Grad.  1 He  uses  some  moments  of  low  order  in  addition  to 
the  usual  ones  representing  p,  q1  and  T.  The  procedure  involves  a gain  in  sim- 
plicity over  the  E-C  theory  but  is  still  quite  complicated.  It  is  basically  a high- 
density  theory  and  is  capable  of  dealing  only  with  boundary  conditions  which  can 
be  specified  in  terms  of  the  moments  appearing  in  the  theory.  The  limitations 
of  any  theory  based  on  the  use  of  a finite  number  of  moments  will  be  brought  out 
in  the  discussions  of  this  series  of  papers. 

The  opposite  limiting  case  (T«  T or  L,  <<  L,c)  has  been  studied  exten- 
sively by  mean  free-path  methods.  Jaffe  ^ has  shown  that  this  case  can  be  treated 

2)  K.  F.  Herzfeld,  Ann.  Physik  23,  465  (1935);  cf.  also  H.  Primakoff,  J.  Acoust. 

Soc.  Am.  14,  14  ^1942);  H.  Tsien  and  R.  Schamberg,  ibid.  18,  334  (1946); 

C.  S.  W.  Chang  and  G,  Uhlenbeck,  Dept.  Eng.  Research,  Univ.  Michigan, 

APL/JHUCM-443,  UMH-3-P,  Feb.,  1948;  M.  Kohler,  Abh.  Braunschweig. 

Wiss.  Gesell.  2,  104  (1950). 

3)  H.  Grad,  Comm.  Pure  Applied  Math.  2,  331  (1949). 

4)  H.  Jaffe,  Ann.  Physik  6,  195  (1930). 
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from  the  point  of  view  of  the  Boltzmann  equation  by  expanding  the  distribution 
function  in  a series  of  inverse  powers  of  the  mean  free  path;  the  first  approxi- 
mation consists  in  neglecting  collisions  completely.  In  contrast  to  the  high 
density  limit,  the  problem  can  be  solved  subject  to  correct  microscopic  bounda- 
ry conditions.  The  method  becomes  extremely  complicated  when  one  attempts 
to  carry  out  higher  approximations  to  give  results  which  go  much  beyond  those 
obtainable  by  simple  mean  free-path  arguments.  For  this  reason  this  method 
has  not  been  applied  to  any  great  extent. 

The  low-density  case  is  of  particular  interest  for  ionized  gases  where 
the  coupling  of  the  material  motions  with  the  electromagnetic  field  plays  an 
essential  and  determining  role.  Ionized  gases  occur  widely  in  nature  and  in  the 

laboratory  with  densities  and  temperatures  covering  a large  range.  In  laborato- 

-3 

ry  discharge  tubes  and  similar  devices  pressures  ranging  from  10  mm  up  to 

4 -2 

10  mm  are  commonplace.  In  the  negative  glow  of  a discharge  tube  at  10  mm 

the  material  is  at  least  singly  ionized;  the  electron  density  is  about  10^/cc  and 

the  mean  free  path  is  of  the  order  of  a few  cm.  At  pressures  of  the  order  of  an 

atmosphere  the  ionization  is  only  partial;  the  electron  density  is  of  the  order  of 

15  y -5 

10  / cc  and  mean  free  paths  are  of  the  order  of  10  cm.  In  nature  the  range 

covered  is  considerably  wider?)  The  interstellar  medium  has  a density  varying 

from  10  ^ gm  /cc  up  to  10  ^ gm/cc.  The  gas  may  be  neutral  as  in  the  low- 

temperature  HI  regions,  or  completely  ionized  as  in  HII  regions  near  very 

luminous  stars.  In  the  outer  corona  of  the  sun  the  degree  of  ionization  is  high 

and  the  electron  density  is  of  the  order  of  10  /cc.  As  we  proceed  in  towards  the 

center  of  the  star,  the  degree  of  ionization  varies  considerably.  In  the  photo- 

spheric  layers  (mainly  responsible  for  the  sun's  continuous  spectrum),  the 

particle  density  is  of  the  order  of  10^/cc  anc*  ionization  is  small.  Near  the 

center  of  the  sun  we  of  course  have  very  much  greater  densities  and  temperatures 


5)  J.  A.  Hynek,  "Astrophysics,  " McGraw-Hill  Book  Co.  , New  York,  1951. 


- 6 - 


and  a very  high  degree  of  ionization. 

We  thus  see  that  it  would  be  very  desirable  to  have  a theory  capable  of 
treating  the  whole  range  from  low  to  high  densities  for  mixtures  of  neutral  and 
ionized  gases.  The  treatment  of  an  ionized  gas  involves  certain  new  elements 
in  addition  to  those  for  neutral  gases;  this  will  be  discussed  in  detail  in  Sec.  2. 
The  region  of  intermediate  density  (T  — T or  L — L ) has  thus  far  been  practi- 
cally unexplored  on  account  of  excessive  mathematical  complications.  The 
procedures  discussed  above  for  the  two  limiting  cases  of  high  and  low  density 
are  not  applicable  in  this  region.  For  similar  reasons,  the  solution  of  definite 
initial  and  boundary  value  problems  with  the  Boltzmann  equation  has  also  not 
been  studied  to  any  extent.  In  order  to  make  some  progress  in  the  investigation 
of  such  problems  it  is  necessary  to  simplify  the  mathematical  procedures  quite 
considerably. 

The  aim  of  this  series  of  papers  is  to  propose  a treatment  of  collision 
processes  in  gases  which  leads  to  just  such  a simple  mathematical  formalism. 
The  main  difficulty  in  handling  the  full  Boltzmann  equation  arises  from  the  com- 
plicated nature  of  the  collision  terms.  In  contrast  to  other  procedures,  ours 
consists  in  replacing  this  troublesome  collision  integral  by  mathematically 
simpler  terms;  these,  however,  are  chosen  so  as  to  conform  to  the  conservation 
laws  for  mass,  momentum  and  energy  and  at  the  same  time  to  represent  certain 
essential  features  of  collisions,  e.g.,  persistence  of  velocity.  This  approach 
makes  possible  a survey  of  the  whole  range  from  low  density  to  high  density 
and  including  the  intermediate  region.  It  leads  to  the  correct  asymptotic  be- 
havior in  the  two  limiting  cases  mentioned  above.  The  mathematical  simplifi- 
cation introduced  with  the  model  enables  one  to  solve  problems  which  are  physi- 
cally more  complex  than  those  soluble  with  the  standard  Boltzmann  equation; 
in  addition,  one  is  able  to  treat  definite  initial  and  boundary-value  problems. 
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Each  collision  term  in  a Boltzmann  equation  consists  of  two  parts.  The 
one  part,  ^ - JJ  fj(^,  t)  *2^2'  ^12  ^ represents  particles  re- 

moved or  absorbed  from  a definite  velocity  range  by  collisions;  the  other  part 
represents  the  particles  emitted  into  that  range  as  a result  of  collisions.  Our 
absorption  term,  -n(x,  t)  k(\T,  x,  t)  f(^,  x,  t),  where  n(x>,  t)  = x",  t)dv 

is  substantially  the  same  as  that  in  the  Boltzmann  equation.  The  emission  term 
(which  is  the  real  source  of  difficulty)  is  replaced  in  our  crudest  model  by  a 
term  representing  a Maxwellian  distribution  of  the  emitted  particles  with  a den- 
sity, mass  velocity  and  temperature  satisfying  mass,  momentum  and  energy 
conservation  requirements.  This  bears  some  resemblance  to  the  "local  thermo- 
dynamic equilibrium"  model  used  in  discussing  the  formation  of  the  continuous 
spectrum  in  stellar  atmospheres.  The  similarity  of  the  Boltzmann  equation  to 
that  governing  propagation  of  radiation  has  been  discussed  by  Jaffe.4^  The 
method  proposed  here  is  quite  different,  however,  and  is  capable  of  considerable 
extension  and  generalization  to  include  the  study  of  physically  complex  problems. 
We  shall  show  in  a later  paper  that  it  may  be  used  to  provide  a method  of  so- 
lution of  the  standard  Boltzmann  equation. 

The  plan  of  this  series  of  papers  is  the  following.  In  the  first  paper  we 
study  systems  which  are,  or  may  be  reduced  to,  one-component  systems;  to 
illustrate  the  technique  we  avoid  mathematical  complications  by  making  the  sim- 
plified assumption  of  constant  collision  time.  This  simple  model  is  applied  to 
the  study  of  dispersion  and  absorption  of  sound  in  a simple  unionized  gas  and 
also  to  the  oscillations  of  an  ionized  gas  with  the  additional  simplifying  assump- 
tion that  the  positive  ions  may  be  replaced  by  a uniform  constant  distribution 
of  positive  charge.  The  physically  most  significant  feature  here  is  the  study  of 
the  transition  from  low-pressure  plasma  waves  to  high-pressure  sound-type 


waves. 
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In  the  second  paper  we  extend  this  treatment  to  a consideration  of  two- 
and  three-component  systems.  This  includes  a discussion  of  the  translational 
dispersion  effect  for  sound  waves  in  unionized  gas  mixtures.  The  further 
effect  on  small  oscillations  arising  from  the  coupling  between  positive  ions 
and  neutrals  will  be  studied  over  the  whole  range  from  low  to  high  pressures; 
the  two  types  of  wave  which  occur  at  low  pressures  are  found  to  merge  into  a 
single  type  of  wave  as  the  pressure  increases  and  the  collision  mechanism  begins 
to  predominate. 

Paper  III  will  deal  with  the  oscillations  of  ionized  gases  in  static  external 
magnetic  fields.  The  range  of  validity  of  the  conventional  magnetohydrodynamic 
theory  will  be  discussed  there.  An  important  point  dealt  with  in  that  paper  is 
the  way  in  which  gaps  in  the  frequency  spectrum,  present  at  low  densities,  dis- 
appear as  the  density  is  increased. 

The  discussions  of  the  first  three  papers  are  confined  to  initial  value 
problems  in  the  theory  of  small  oscillations.  The  accent  there  is  substantially 
on  the  study  of  the  mechanism  of  co-operation  in  ionized  and  unionized  gases. 

In  Paper  IV  we  abandon  the  assumption  of  constant  collision  time  and  examine 
physically  more  realistic  models.  The  connection  of  our  kinetic  equations  with 
the  standard  Boltzmann  equation  is  discussed  in  more  detail.  In  subsequent 
papers  we  shall  study  particular  boundary-value  problems  and  some  nonlinear 
processes. 

2.  EARLIER  TREATMENTS  OF  SYSTEMS  OF  CHARGED  PARTICLES 

The  theory  of  oscillations  in  ionized  gases  is  still  in  a rather  tentative 
state.  However,  a number  of  features  can  be  understood  at  least  qualitatively. 

An  essential  difference  from  the  theory  for  neutral  gases  is  that  approximately 
undamped  waves  can  be  propagated  not  only  at  high  pressures  (t  » Tc) 
but  also  at  low  pressures  (T  « Tc). 
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The  mechanisms  responsible  for  the  co-operative  behavior  of  the  medium 
are  different  in  the  extreme  cases  of  low  density  and  high  density.  ^ At  low 
densities  the  collisions  are  only  of  secondary  importance  and  the  forces  acting 
to  change  the  state  of  motion  of  a particle  are  electromagnetic  in  character. 
Energy  and  momentum  are  transferred  from  one  group  of  particles  to  another 
through  the  intermediary  of  the  electromagnetic  field  in  a way  now  familiar 
from  the  study  of  microwave  devices.  It  is  important  to  note,  however,  that 
the  transfer  of  energy  by  the  high-frequency  components  of  a pulse  is  a very 
slow  process  in  low-pressure  plasmas.  If  one  neglects  the  effect  of  random 
thermal  motion,  the  group  velocity  of  the  waves  is  zero  and  a disturbance 


would  remain  completely  localized.  Including  the  effects  of  thermal  motions, 

the  group  velocity  of  high-frequency  waves  is  approximately  v = (3pa/&»  )a, 

/ 8 P 

where  a =(kT/m)  ' ; this  is  much  smaller  than  the  sound  velocity,  5a/3,  for 

waves  whose  wave  number  p is  smaller  than  the  inverse  of  the  Debye  length 

(a/o>p).  In  addition,  waves  with  p££  a/uip  are  very  heavily  damped. 

Let  us  now  turn  attention  to  attempts  which  have  been  made  to  treat  the 

oscillations  of  ionized  gases  quantitatively  by  the  methods  of  kinetic  theory. 

It  is  clear  that  the  long  range  of  the  Coulomb  force  implies  that  the  concept  of 

a definite  binary  collision  will  be  inadequate  except  for  very  close  encounters. 

Instead,  every  charged  particle  is  continually  interacting  with  all  the  particles 

of  the  assembly.  An  important  theoretical  step,  making  possible  the  approxi- 

7) 

mate  treatment  of  many  properties  of  ionized  gases,  was  taken  by  Vlasov.  ' 

8) 

Langmuir  ' had  already  pointed  out  that  the  Debye  length  is  characteristic  of 
the  dynamic  behavior  as  well  as  of  the  static  properties  of  an  ionized  gas. 


6)  D.  Bohm  and  E.  P.  Gross,  Phys.  Rev.  75,  1851,  1864  (1949). 

7)  A.  Vlasov,  J.  Phys.  U.S.S.R.  9,  25,  130  (1945). 

8)  I.  Langmuir,  Proc.  Nat.  Acad.  Sci.  U.S.  14,  627  (1928). 
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Processes  involving  characteristic  lengths  longer  than  the  Debye  length  call 
into  play  the  co-operative  behavior  of  the  particles  of  an  ionized  gas.  Vlasov 
assumed  that  apart  from  close  collisions,  a particle  moves  in  an  average 
electric  field  arising  from  the  other  particles  and  proposed  to  describe  the 
system  by  a one -particle  distribution  function  f(v,  x,  t).  He  wrote  for  the 
description  of  longitudinal  disturbances, 

7E  + (v.v)f._  ^=^1  , (1) 


V'  E = {n rJi  d?} 


where  nQ  is  the  density  of  the  uniform  continuous  positive  background.  The 
force  eE  arising  from  the  charge  concentrations  of  the  particles  themselves 
has  been  inserted  in  the  position  in  the  kinetic  equation  generally  reserved  for 
body  forces.  Vlasov  termed  this  the  integral  interaction  term;  it  is  designed 
to  treat  the  forces  on  a given  particle  arising  from  other  particles  at  distances 
greater  than  a Debye  length.  The  term  | represents  the  instantaneous 
change  of  the  distribution  function  arising  from  closer  collisions.  The  question 


then  arises  as  to  the  explicit  form  to  take  for 


The  very  close  encounters 


at  a distance  less  than  the  interparticle  distance  can  be  treated  by  the  binary 
collision  term  of  Boltzmann  type.  For  the  more  distant  collisions  such  an  ap- 
proach is  less  certain.  It  is  important  to  note  that  in  a low-pressure  plasma 
the  Debye  length  is,  in  general,  many  times  the  interparticle  distance.  Thus, 

in  a collision  between  two  particles  there  are,  in  general,  many  other  particles 

9) 

between  them.  In  spite  of  this,  many  authors  ' have  treated  various  processes 
in  ionized  gases  using  the  binary  collision  approach  even  for  these  distant 


9)  D.  Gabor,  Z.  Physik  84,  474  (1933);  L.  Landau,  Physik.  Z.  Sowjetunion 
1£,  154  (1936);  D.  B ohm  and  L.  H.  Aller,  Astrophys.  J.  105,  131  (1947); 

R.  Cohen,  L.  Spitzer,  Jr.,  and  P.  McR.  Routly,  Phys.  Rev.  80,  230  (1950). 
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collisions.  They  treat  the  distant  collisions  as  involving  small  momentum 
transfers  and  expand  the  distribution  functions  in  the  collision  integral  term 
in  powers  of  the  momentum  transfers.  Terms  very  similar  to  those  given 
by  the  Fokker -Planck  equation  are  obtained.  This  procedure  is,  of  course, 
approximate  since  it  attempts  to  replace  many  body  collisions  by  a series  of 
binary  collisions.  The  equation  neglects  the  "Onsager  correction”  and  other 
correlation  effects;  these  require  a more  refined  description  in  terms  of  many 
particle  distribution  functions.  However,  in  the  present  primitive  state  of  the 
theory  it  is  worthwhile  exploring  the  consequences  of  this  approach.  In  a later 
paper  we  shall  discuss  the  application  of  this  method  to  the  plasma  oscillation 
problem  in  some  detail.  The  Fokker-Planck  approach  has  already  been  used 
in  the  study  of  plasma  oscillations  by  Logunov.  ^ He  obtains  a damping  and 
dispersion,  but  his  method^  is  incapable  of  treating  high-pressure  plasmas  be- 
cause of  an  inadequate  representation  of  the  collision  terms. 

3.  KINETIC  MODELS  OF  COLLISION  PROCESSES 

The  equations  described  above  for  neutral  and  ionized  gases  are  compli- 
cated to  handle  in  practice  because  of  the  intractable  nature  of  the  Boltzmann 
binary  collision  term.  We  shall  now  discuss  some  simple  kinetic  models  which 
permit  of  exact  mathematical  treatment  including  the  solution  of  definite  boundary- 
value  problems.  A more  detailed  discussion  of  the  validity  of  our  representation 
of  collisions  and  its  relation  to  the  standard  Boltzmann  equation  will  be  given  in 
a subsequent  paper.  In  this  section  of  the  present  paper  we  introduce  the  simplest 
kinetic  model;  we  give  the  underlying  motivation  for  the  form  of  this  model  from 
one  point  of  view. 

In  many  other  kinetic  problems,  for  example,  the  effect  of  electron  col- 
lisions on  the  propagation  of  radio  waves  in  the  ionosphere,  it  is  convenient  to 


10)  A.  A.  Logunov,  Zhur.  Eksptl.  Teoret.  Fiz.  2£,  458  (1950). 
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avoid  the  complexities  of  the  Boltzmann  equation  by  using  a mean  free-path  treat- 
ment. One  replaces  the  collision  integral  by  a relaxation  term  of  the  form 


6f 

Tt 


fo  - f(  V,  X,  t) 

r l v) 


(3) 


where  T(  v)  is  a velocity  dependent  collision  time.  This  expresses  the  fact 
that  collisions  tend  to  relax  the  distribution  function  to  an  equilibrium  value  fQ. 
We  illustrate  our  discussion  of  collision  models  by  referring  to  oscillatory 
problems  where  a characteristic  time  enters  in  a natural  way.  The  relaxation 
model  then  describes  the  destruction  of  phase  of  an  ordered  motion  on  collision 
and  leads  to  a damping  frequency  of  order  l/T  in  the  amplitude,  where  T is 
some  suitable  average  collision  time.  This  type  of  model  has  the  defect  ^ that 
charge  is  not  conserved  instantaneously  but  only  over  a cycle.  It  is,  however, 
easy  to  remedy  this  at  least  in  the  case  of  constant  collision  time  by  taking  for 
the  collision  term 


= . 1+  £lLL±L  f (S) 

6t  T T o^v; 


(4) 


where  n(x,  t ) = J'  f dv  is  the  fluctuating  density.  Thus  particles  in  a range  dv 
about  velocity  v are  absorbed  at  a rate  proportional  to  the  number  f ( v,  x,  t)  at 

(x,  t)  and  reemitted  at  a rate  proportional  to  the  density  at  (x,  t)  and  with  a 

* 

Maxwellian  distribution  of  velocities.  We  may  regard  this  as  a model  represent- 
ing electron  neutral  collisions.  (This  representation  will  be  studied  in  more 
detail  in  Paper  2,  where  we  treat  the  coupling  of  charged  and  uncharged  systems 

more  adequately.)  Upon  integrating  over  velocities  we  find  J -gl  dv  = 0 so 

c 

that  the  number  of  particles  is  instantaneously  conserved  in  collisions  in  contrast 
to  the  earlier  case  where  the  conservation  holds  only  when  averaged  over  finite 
times.  There  is,  however,  no  instantaneous  conservation  of  momentum  or 
11)  E.  P.  Gross,  Phys.  Rev.  82,  232  (1951). 

* For  a discussion  of  this  model  from  a different  point  of  view,  Ref.  6,  p.  1864. 
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energy.  In  fact,  we  have 


/*  -+  6f  I , -*  m f -*  c /‘mv^  6f  I , -*■ 

Jmvw\  dv  = - — J vf  dv,  J -j-  w I dv  = 

•c  c 

This  represents  a gain  or  loss  of  total  energy  and 


momentum  during  collision,  depending  on  the  phase  of  f. 

Thus  far  we  have  discussed  a model  which  can  represent  only  electron- 
neutral  or  ion-neutral  collisions  in  which  the  collision  term  is  linear  in  the 
electronic  distribution  function.  The  linear  feature  would  also  be  present  in  a 
more  exact  study  of  the  collisions.  However,  in  considering  electron-electron 
or  ion-ion  collisions  which  are  important  in  the  high-pressure  ionized  gas,  we 
will  expect  to  find  collision  terms  quadratic  in  the  electron  and  the  ion  distribution 
functions.  (In  the  study  of  small  amplitude  disturbances  such  terms  may  be 
linearized,  but  we  now  seek  a general  formulation.)  We  shall  construct  a simple 
model,  quadratic  in  f which  reproduces  essential  features  of  the  collisions  with- 
out introducing  the  troublesome  Boltzmann  collision  integral.  Let  us  write 


df 


2 y 4 v 

■.  t)  +2_ti_y  f 


(5) 


2 

_ m v 

F = ( , ™ )3/2  e"  *kt° 

o 

Thus  J F dv  = 1.  We  have  taken  the  rate  at  which  particles  leave  a particular 
velocity  range  at  (x,  t)  as  proportional  to  f(v,  x,  t)  and  to  the  density  n(x,  t). 
The  number  of  particles  emitted  at  ( x,  t)  equals  the  number  absorbed  at  (x,  t); 
this  number  is  proportional  to  the  square  of  the  density  and  is  emitted  with  a 
Maxwellian  velocity  distribution.  H = nQ/< r is  a collision  frequency,  here 
assumed  independent  of  velocity.  It  is  easily  verified  that  the  total  number  of 
particles  is  instantaneously  conserved  by  the  collision  processes,  while  momen- 
tum and  energy  are  not. 

It  is  now  worthwhile  to  examine  the  dispersion  relation  for  small -ampli- 
tude oscillatory  processes  governed  by  Eq.  (5).  For  brevity,  we  adopt  the 
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approach  of  looking  for  plane-wave  solutions  of  the  type  exp  (i  (px  - «t) ) . We 
shall  show  later  that  this  procedure  is  valid  for  the  treatment  of  certain  types 
of  initial  value  problems  for  finite  collision  time  and  for  waves  of  large  wave- 
length. We  write 


f = fQ  (1  + ♦) 


f = n F 
o o 


(7) 

n (x,  t)  = nQ^(l  + V ) , (8) 

where  and  V are  dimensionless  quantities  small  compared  to  unity.  We  as- 
sume 4*.  V , and  E oscillate  as  exp  (i  (px  - ut)  ) and  neglect  products  of  these 
variables.  Then  Eq.  (5)  becomes 


i (pu  - <->)  <j>  + pj,  u = X.  ( V - 4>) 
o 

V = J F <j>  dv 
4iren 


We  also  have 


E = - 


IP 


Eliminating  E we  find  the  distribution  function 


^ V.  + i (pu  - w)  j 

and  the  dispersion  relation 

1 = 


X.  + 


4im  e 
o 


TT~  Ip 


f F dv 
J V.  + i (pu  - u) 


X + 


4irn  e 
o 

TT~ 


u 

ip 


(9) 

(10) 
(H) 

(12) 

(13) 


If  one  expands  in  powers  of  ipu/(X  - iw)  and  restricts  oneself  to  the  study  of  long 


wavelength  one  finds 


2 2 
w = u> 

P 


kT 


iu»\  - 


o 2 
- P 


X.  - iu> 


(14) 


(X  - ic*>) 

We  conclude  from  this  relation  that  the  oscillations  are  strongly  damped  at  high 

pressures.  For  the  case  of  a charged  gas  the  complex  frequency  is  approximate 
iX  / 7 \ 1/2 

ly  w zz  - Ar-  +1  - -r  + J • this  is  of  the  order  of  the  collision  frequency  for  X 
2 * ' 4 p / ' .2 


small  and  always  corresponds  to  damping.  For  a neutral  gas,  w = ( 

2 kT  2 V P 

and  so  w — - iwX  + ——  p -r — 

m X - iu 


4trn  eL  . / , 
o \1/  2 


m 


) 


; thus  in  the  limit  of  very  high  collision  fre- 


quency one  obtains  highly  damped  oscillations  instead  of  undamped  sound  waves. 
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We  thus  see  that  this  model  is  incapable  of  yielding  correct  results  at 
the  high-pressure  limit?  What  is  the  origin  of  this  failure  . As  discussed 
earlier  and  in  Ref.  6,  the  processes  of  organization  are  very  different  for  a 
plasma  and  for  a sound  wave.  In  particular,  one  must  recognize  the  velocity 
organization  characteristic  of  a sound  wave.  Thus  a particle  with  excess 
velocity  at  a point  makes  collisions  which  slow  it  down  to  the  flow  velocity  at 
the  point  in  question.  Similarly  particles  which  are  moving  too  slowly  tend  to 
have  their  speed  increased  to  that  of  the  flow  velocity.  The  collision  models 
studied  thus  far  fail  to  incorporate  this  feature  so  that  it  is  not  surprising  that 
they  do  not  yield  sound-type  oscillations  in  the  high-density  limit. 

We  now  propose  a simple  model  which  succeeds  in  giving  an  adequate 
description  of  both  the  high-  and  low-density  plasma.  For  a one-component 


system,  one  writes 


• v> « - 1- u- '(*■ *.  t) 


<&  = ( 


n \ 3/2  I m # -*■  -*• 

— T — ) exp — ( v - q 

2irkT  ( x,  t ) L 2kT  ( x,  t ) 

J <J>  d v = 1 

n ( x,  t ) = J f ( v,  x,  t ) d v 

q (x,  t)  J v f ( v,  x,  t)  dv 

n(£.  tr 


(15) 

(16) 

(17) 

(18) 


= 1 

m n J ' n ' 

V-  E = 41re|no-/fd7j. 

q(x,  t)  and  T (x,  t)  define  the  flow  velocity  and  temperature  at  x and  t. 


(19) 

(20) 

We  have 


assumed  the  reemitted  particles  at  ^ x,  t)  to  emerge  with  a Maxwellian  distri- 
bution centered  about  the  flow  velocity  q and  for  a temperature  T ( x,  t).  It  is 


easily  verified  that,  with  the  definitions  of  q and  T given,  the  collision  terms 
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instantaneously  conserve  particle  number,  momentum,  and  energy.  The  ki- 
netic equation  (15)  is  a nonlinear  integro-partial  differential  equation.  How- 
ever, it  is  considerably  simpler  than  the  standard  Boltzmann  equation  since 
the  distribution  function  enters  into  the  collision  terms  of  Eq.  (15)  in  a simple 
way;  as  f,  Ji  d v,  J v f dv  and^ v^  f dv.  Thus  one  finds  a solution  of  Eq.  (15) 
in  terms  of  the  undetermined  functions  q(x,  t),  n(x,  t),  T(  x,  t),  E(  x,  t),  and 
then  inserts  this  solution  into  Eqs.  (17)  to  (20).  This  establishes  the  con- 
ditions of  compatibility  which  serve  to  determine  the  solution  completely.  For 
the  particular  case  of  small  amplitude  oscillations  the  procedure  yields  a dis- 
persion relation  specifying  the  connection  between  frequency  and  wavelength  in 
its  dependence  on  collision  frequency,  plasma  frequency,  temperature,  etc. 

The  use  of  conditions  of  compatibility  bears  a certain  resemblance  to  the 
first  step  in  the  Enskog -Chapman  method  of  solving  the  Boltzmann  equation. 
However,  the  methods  are  quite  different.  In  the  E-C  technique  one  assumes  a 
locally  Maxwellian  form  for  the  distribution  function  in  first  approximation  and 
then  determines  the  density,  velocity  and  temperature  entering  in  this  choice 
by  conditions  of  compatibility;  here  we  change  the  collision  term  itself  and 
solve  for  the  distribution  function  rigorously.  The  solution  will,  in  general,  not 
have  the  form  assumed  in  the  Enskog -Chapman  method.  Our  approach  has  the 
virtue  that  one  can  solve  definite  initial  and  boundary -value  problems  wholly 
within  the  framework  of  a microscopic  formalism.  The  model  we  have  described 
involves  the  assumption  of  a collision  time  independent  of  velocity.  In  Paper  4 
we  shall  show  that  this  is  not  an  essential  restriction  and  shall  discuss  the  gener- 
alizations to  take  into  account  velocity-dependent  collision  times  as  well  as 
others;  in  the  present  paper  we  shall  follow  the  consequences  of  Eqs.  (15)  to  (20) 
inasmuch  as  it  appears  to  be  the  simplest  model  capable  of  yielding  physically 
satisfactory  results  throughout  the  entire  range  of  pressures. 
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There  exists  a number  of  studies  of  the  oscillations  of  ionized  gases 

121 

using  transport  equations.  For  example,the  Thomsons  ' and  subsequently 
13) 

Bailey  'base  their  discussions  on  Maxwell's  equations  of  transfer.  As  ap- 
plied to  a one -component  system  of  charged  particles  these  equations  are  those 

of  continuity,  of  momentum  transfer  and  Poisson's  equation.  When  no  static 

2 2 kT  2 

magnetic  field  is  present  one  obtains  the  dispersion  relation  i*>  = + (-^-  )p  • 

It  is  not  clear  from  these  transport  treatments  whether  the  results  are  meant 

to  apply  at  low  pressures  or  high  pressures;  in  fact,  they  could  be  valid  only  at 

high  pressures.  Even  there,  however,  they  are  inexact,  as  they  correspond  to 

an  isothermal  treatment.  Thus  , in  terms  of  the  model  discussed  above,  the 

temperature  T is  treated  as  a constant  and  energy  balance  is  not  maintained 

instantaneously  by  collisions;  we  shall  obtain  Thomson's  results  from  our 

model  at  high  pressures  in  the  isothermal  approximation.  The  high-pressure 

nature  of  the  Thomson  approach  for  a one-component  system  is  brought  out  in 

14) 

an  alternative  treatment  by  Linder  ' who  adds  a pressure  term  to  the  equation 
of  motion  of  a particle  and  with  an  isothermal  assumption  obtains  the  same 


result  as  Thomson.  At  low  pressures  the  correct  dispersion  relation  as  found 
by  an  exact  treatment  of  the  kinetic  equation  is  u = — p , differing 

from  that  at  high  pressures  by  the  factor  3.  We  expect  then  that,  in  an  iso- 
thermal treatment,  there  will  be  a change  in  the  frequency  corresponding  to  a 
given  wavelength  as  one  goes  from  low  to  high  pressures.  This  is  indeed  the 
case  as  will  be  shown  with  our  collision  model;  in  addition,  we  find  an  accom- 
panying absorption  which  reaches  a maximum  value  when  the  collision  frequen- 
cy is  of  the  order  of  magnitude  of  the  plasma  frequency. 

12)  J.  J.  and  G.  P.  Thomson,  "Conduction  of  Electricity  through  Gases,  " 
Cambridge  University  Press,  London,  1933,  Vol.  2,  p.  353. 

13)  V.  A.  Bailey,  Phys.  Rev.  78,  428  (1950). 

14)  E.  G.  Linder,  Phys.  Rev.  49j  753  (1936). 
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We  have  thus  far  confined  discussion  of  the  transport  method  to  iso- 
thermal models.  It  is  possible  to  generalize  the  transport  procedure  to  in- 
clude temperature  variations  by  applying  a moment  method  to  the  Boltzmann 
equation.  ^ One  finds  a dispersion  relation  + y p^  kTQ/m.  For  the 

uncharged  case  this  gives  the  correct  "adiabatic"  law  of  sound-wave  propa- 
gation u = (5/3)^^  (kTQ/m)^^p.  We  shall  see  that  the  result  is  also  correct 
for  an  ionized  gas;  the  transport  approach  owes  its  validity  at  high  pressures 
to  the  fact  that  momentum  and  energy  are  collisional  invariants.  Existing 
transport  treatments  are  certainly  incorrect  quantitatively  in  the  low-pressure 
region  and  even  fail  to  give  a qualitatively  adequate  description  when  magnetic 
fields  are  present.  ^ For  the  treatment  of  the  effects  of  heat  conduction  and 
viscosity  which  depend  on  the  mean  free  path  a kinetic  approach  of  the  type 
here  presented  is  needed. 


4.  ISOTHERMAL  APPROXIMATION 
The  simplest  collision  model  yielding  both  low-pressure  plasma  waves 
and  also  sound  waves  at  high  pressures  may  be  obtained  by  taking  T constant  in 
Eqs.  (15)  to  (20).  This  is  not  entirely  realistic  since  it  yields  the  Newtonian 
isothermal  sound  speed  at  high  pressures.  But  the  fact  that  the  number  of  vari- 
ables is  one  less,  simplifies  the  calculations  considerably.  For  the  sake  of 
ease  of  exposition  we  shall  discuss  this  case  in  detail  and  only  sketch  the  treat- 
ment of  the  more  general  case  including  temperature  variations.  The  isothermal 
case  is  also  of  some  interest  in  itself  since  the  work  of  Thomas  and  Bailey  is 
based  on  the  isothermal  approximation.  The  results  of  main  physical  interest 
are  presented  in  Figs.  1 to  5.  Figures  1 and  2 show  the  dispersion  and  absorption 
of  a one-component  charged  system  as  a function  of  collision  time.  Figure  5 shows 
the  absorption  and  the  translational  dispersion,  effect  in  a neutral  monatomic  gas. 
The  equations  describing  the  system  of  particles  in  the  approximation  of 
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this  section  are  Eqs.  (15)  to  (18),  (20).  We  do  not  use  Eq.  (19)  but  set  the 
temperature  equal  to  a constant  Tq,  the  equilibrium  temperature  of  the  gas. 

As  we  are  concerned  with  small  amplitude  oscillations  we  shall  neglect  quanti- 
ties quadratic  in  the  oscillating  variables  and  arrive  at  a set  of  linear  equations. 
The  consistency  of  the  linearization  procedure  can  be  checked  by  comparing 
the  terms  neglected  (as  estimated  with  the  linear  approximation)  with  the  terms 
retained.  For  the  case  of  low-pressure  plasma  oscillations  a breakdown  of  the 
linear  approximation  is  closely  connected  with  the  presence  of  a large  number 
of  trapped  particles  moving  with  a velocity  close  to  the  wave  velocity.  These 
particles  play  a role  quite  different  from  other  particles  and  make  possible 
types  of  solutions  outside  the  scope  of  a linear  approximation. 


f ( v,  x,  t ) = nQ  F { 1 + $ | (21) 

n (x,  t)  = nQ  | 1 + yj.  . (22) 

Then  <J>  and  V are  dimensionless  quantities.  Equations  (15)  to  (20)  then  become 

Z/(x,  t)  = /«j.  F dv  (23) 

q(x,  t)  =y  v < j>  F d v (24) 

^ ’ E = -4ire  n qV  . (25) 

In  the  linearized  theory,  <}»,  V , q,  E are  quantities  of  first  order  of  smallness. 
We  consistently  neglect  all  products  and  higher  powers  of  these  quantities.  On 
expanding  the  exponential  in  Eq.  (16)  (neglecting  q^  and  higher  powers),  we  find 


<J>  = F / 1 + 


The  linearized  kinetic  equation  then  becomes 


!£+(v.  y)<t>  + ^-?  = \(i/  -♦+1?J2-v-  I)  (27) 

° o o 

Equations  (21)  to  (27)  form  the  basis  for  the  further  analysis  of  this  section. 

15) 

Landau  ' has  emphasized  the  lack  of  rigor  in  the  procedure  of  looking  for 


15)  L.  D.  Landau,  J.  Phys.  U.S.S.R.  10,  25  (1946). 
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solutions  of  Eq.  (27)  of  the  type  where  all  quantities  have  a time  and  space 
variation  e*  One  must  consider  also  the  solutions  of  the  homogeneous 

parts  of  Eqs.  (20)  and  ^27)  and  take  the  correct  linear  combination  in  order  to 
solve  a definite  initial  or  boundary -value  problem.  In  this  paper  we  examine  a 
one-dimensional  initial  value  problem,  the  propagation  of  a pulse.  More  general 
boundary- value  problems  will  be  treated  in  subsequent  papers.  Suppose  that  at 
times  t < 0 the  particle  density  is  uniformly  nQ  throughout  unbounded  space 
except  for  a small  region  between  - xq  and  xq  confined  between  two  walls,  where 
the  density  is  n + v . At  time  t = 0 the  walls  are  removed  and  we  ask  about 
the  subsequent  behavior  of  the  system.  To  solve  the  problem  the  distribution 
function  at  t = 0 must  be  specified.  This  initial  distribution  depends  on  the  par- 
ticular (macroscopic)  mode  of  preparation  of  the  system.  If,  for  example,  the 
present  system  has  been  prepared  by  compressing  the  gas  in  ( - xq,  +xq)  slow- 
ly over  a time  very  long  compared  to  the  mean  time  between  collisions,  the 
initial  distribution  function  will  be  very  nearly  Maxwellian.  For  the  case  of 
sound  waves  it  is  clear  what  will  happen.  Two  pulses  are  propagated  with  sound 
speed  in  opposite  directions  parallel  to  the  x axis;  the  shape  of  each  pulse  is 
maintained  except  for  the  small  absorption  and  dispersion  effects.  In  a low- 
pressure,  ionized  gas  the  pulse  will  move  slowly  but  will  alter  its  shape  appre- 
ciably because  of  the  high  dispersion  . 

We  take  as  initial  conditions  at  t = 0 

<t>  ( v,  x,  0)  = VQ  F ( v ) for  - xq  ^ x < xq 
= 0 otherwise 

, ) tZ8) 


V - v_ 


for  - x < x < x 
o o 


J 


- o otherwise 

The  mass  velocity  qis  zero  throughout  all  space  at  t = 0. 

The  special  problem  is  typical  of  a large  class  of  problems  which  may  be 
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solved  by  first  performing  a Fourier  expansion  in  unbounded  x space.  We  write 


4>  ( v,  x.  t)  =f$ p ( v,  t)  e 1 p ’ x dp 


, etc. 


From  Eqs.  (23)  to  (27)  we  find  the  following  equations  for  the  Fourier  trans- 


forms, 


eE  • v 
P 


jr  + iP‘  — ft-  = ^ UP ' *P  + ft  * ' V ^30) 

♦ O f*  + Q + 

V (t)  =/“<j>  (v,  t)  F dv  (31) 

qp  (t)  - J v 4>p  F dv  (32) 

ip  ■ E = - 4iren  V (t)  (33) 

r £ op'7  . v ' 

These  equations  are  a determined  set  if  one  specifies  cj>p(  v,  t)  at  time  t = 0. 

Thus  we  must  find  the  solutions  of  Eqs.  (30)  to  (33)  for  which  each  Fourier 

component  satisfies  the  appropriate  initial  condition.  One  then  superposes  the 

<|>  to  find  <(>  ( v,  x,  t ) . 

For  example,  in  the  problem  discussed  above 


+x  >oo  + oo 


•'  •/;  /. 


V F e " 1 P * X dx 
o 


sin  p x _ 

= 2 V — 6 (p  ) 6(p  ) F 

op  ' y ' ' 9 


sin  p x 

*'p(°)  =2  Va  p^-2-  6(Py)6(PZ) 


where  6(x)  is  the  Dirac  delta  function.  Let  us  suppose  that  V (t)  has  the  value 

V0)  +i  p 

1/  (t)  = — ^ — ( e 1 + e ).  For  oscillations  of  a one-component  plasma  with- 
out thermal  motions,  w = w is  independent  of  p and  the  disturbance  remains 

P I kTo \ ^ 

localized.  For  isothermal  sound  waves,  w = p(  — — J and  we  obtain  two  oppo- 
sitely directed  pulses.  For  long  wavelengths  and  low  pressures,  the  plasma, 

3 kT  2 

including  thermal  effects,  obeys  the  dispersion  relation  w ci  u + -g  -E- 

P ^ ^r» 
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To  investigate  the  properties  of  a density  pulse  we  form 


r ■ •*  sin  p x 

nit)  - „oy>-*  -Jss. 


6(py)  6(pz) 


iw  t -1 
p 

e r e 


, kT  , 
.3  o 2, 

•K  P t 

2 m r 


+ conj . comp 


.]• 


The  detailed  discussion  of  this  integral  is  complicated  and  will  not  be  given 
here.  In  the  limit  of  long  times,  l/(x,  t)  tends  to  zero  throughout  all  space  in- 
dicating the  spreading  of  the  initially  sharp  pulse  due  to  dispersion.  As  we 
shall  see  later,  in  reality  the  dependence  of  t ) is  more  complicated  than 
that  assumed  here;  for  pulses  with  widths  large  compared  to  the  Debye  length 
the  general  features  are  unaltered. 

We  solve  Eqs.  (30)  to  (33),  following  Landau's  procedure,  by  performing 
a Laplace  transformation  with  respect  to  the  time  variable. 

We  write 

V o- = / t) e_,rt  dt  (35) 

->  >/o  i 


for  the  image  of  <j>^  ( v,  t).  For  the  inverse  transformation  we  have 

»c+ioo 


V 


ft  C+lOO 

J C-lOO 


<rt  , 
e dcr 


(36) 


where  the  integration  is  carried  out  along  a line  in  the  right  half  plane  of  the 
complex  variable  (r.  For  the  case  of  waves  propagated  along  the  x axis  Eqs.  (30) 
to  (33)  become 

i 7ft  (E  ) + X y + \ r-7p — u(q  ) + 4>  ( v,  0) 

kT  ' p,  <rx  p,  <r  kT  xnp,  cr'x  Yp'  ' 


1 (X.  + tr  + ipu)  <(> 


p,  o- 


F dv 


o 

VP>  * = P>  - 

(vA7,v»fd? 

♦ £ 

ip(En  = -4iren  V 

P>  x o p,  or 


(37) 

(38) 

(39) 

(40) 


Here  $^(v,  0)  is  the  Fourier  transform  of  the  initial  velocity  distribution.  By 
multiplying  Eq.  (37)  by  F and  integrating  over  all  velocities  one  finds 


<r  v + ip(q  ) = V (0) 

"p,  tr  *^Mp,  <r'x  "p' 


(41) 
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where  1/^ ( 0 ) = J' $ (v,  0)  F dv. 


Using  this  equation  together  with  Eq.  (40),  we  eliminate  (q  ) and  (E  ) 


p,  cr  x 


from  Eq.  (37)  and  are  left  with  the  set  of  equations 

, 2 


(K  + cr  + ipu)  * =(X-^-X  + I^u)  V +♦  (OJ  + ^Sr 

v o r r o 

"p 

■> 


Solving  for  Z/p>  ^ we  find  u V (0) 

— 4-+y°» 

Vp.  0-  = 2 a . (44) 

"♦  T F d^-  f . Im  ucr  1 

XTTTTpu  | “ FF^  ip  + ip  u J 

2 2 2 

Here  kQ  = «p/a  • The  expression  for  <j>p  ^ is  obtained  by  inserting  this  value  of 

j,  into  Eq.  (42).  Equations  (42)  and  (44)  thus  provide  a complete  solution  to 
p,  cr 

thcTinitial -value  problem  since  we  have  merely  to  invert  the  Laplace  transform 

to  find  the  dependence  of  the  density  and  velocity  distribution  on  time.  However, 

the  integrations  are  along  paths  in  the  right  half  plane  of  the  complex  variable  cr 

and,  aside  from  being  difficult  to  evaluate,  yield  little  physical  insight  into  the 

behavior  of  the  solution.  We  therefore  note  that  Eqs.  (42)  and  (44)  are  defined 

only  for  <r  in  the  right  half  complex  plane  but  that  it  is  possible  to  effect  an 

analytic  continuation  of  the  functions  into  the  left  half  plane  of  a by  modifying  the 

path  of  integration  in  the  integrals  with  respect  to  velocity.  We  integrate  in  the 

complex  u plane  along  a path  C which  passes  below  the  singularity,  X.  + <r  + ipu  = 0 

or  u = - — • The  two  integrals  entering  in  the  numerator  and  denominator 

of  Eqs.  (44)  are  then  entire  functions  in  the  complex  cr  plane. 

With  the  path  of  integration  C,  the  Laplace  transform  of  the  density  y 

P.  cr 

will  have  poles  at  values  of  cr  given  implicitly  by 

“/c  = 0 (45> 


V0-  = 
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There  will  be  an  infinite  number  of  such  poles,  all  lying  in  the  left  half  plane  of 
<r.  If  one  deforms  the  contour  into  the  left  half  plane  in  Eq.  (45),  one  finds  that 
the  density  after  a long  time  t is  given  very  nearly  in  terms  of  the  poles  at  the 
residues  of  V , <r.  For  shorter  times  the  remaining  parts  of  the  contour  give 
a significant  contribution  and  there  is  in  fact  no  advantage  in  deforming  the  con- 


tour in  the  above  way.  For  t large  then 

u 1/(0) 


'(*)-/ 


~ j -►  / Xm  £ 

F dv<FT 


IP 


+ 4>p(  o) ) 


(T.t 
2 e J 


° ^ J 

since  the  numerator  in  Eq.  (44)  is  independent  of  <r.  I/(t)  may  thus  be  repre- 

sented approximately  as  a sum  of  damped  oscillating  terms.  The  asymptotic 
behavior  for  t large  is  determined  by  that  root  of  the  dispersion  relation 
Eq.  (45)  which  has  the  largest  real  part. 

The  dispersion  relation  Eq.  (45)  differs  from  that  obtained  by  the  e^X 
substitution  in  that  the  path  of  integration  is  not  necessarily  along  the  real 
axis.  We  shall  see  shortly,  however,  that  the  substitution  analysis  gives  a 
correct  account  of  the  behavior  at  long  times  when  the  damping  is  not  very  great. 
The  behavior  of  the  velocity  distribution  <j>  ( v,  t)is  different  from  that  of 

the  density  V ( t ) since  the  Laplace  transform  of  d>  has  an  additional  pole  at 

P P,  o-  r 

a value  of  <r  given  by 


X + <ro  + ipu  = 0 . (46) 

In  this  expression  the  velocity  u is,  of  course,  real  since  we  are  investigating 
the  behavior  of  the  distribution  function  for  particles  of  a given  real  velocity  u. 
The  asymptotic  behavior  of  <$>(  v,  t)  depends  essentially  on  the  relative  values 
of  the  real  parts  of  the  pole  given  by  <tq  = - ipu  - X and  of  the  pole  <r ^ with  largest 
real  part  satisfying  the  dispersion  relation  Eq.  (45).  It  is  easy  to  see  that  the 
dispersion  relation  will  be  the  determining  factor  when  the  damping  frequency  is 
smaller  than  the  collision  frequency  X.  For  this  case  both  the  density  and  the 
distribution  function  (for  any  value  of  the  velocity)  have  the  same  time  dependence 
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<r  i 

e .We  shall  call  this  the  "collision-damping"  case.  From  our  later  ex- 
pression Eq.  (68)  for  the  frequency  as  a function  of  wavelength  and  collision 
time,  we  see  that  Re  = w.  = - [^2  + i)J  [pV/<-p]  for  waves  larger  than 
either  the  Debye  length  or  the  mean  free  path.  We  see  that  the  damping  fre- 
quency is  indeed  smaller  than  the  collision  frequency  for  these  waves.  This 
condition  has  the  further  consequence  that  the  contour  C may  be  taken  along  the 
real  axis  so  that  the  e ^pX  " substitution  analysis  will  lead  to  correct  results 
for  the  behavior  at  large  values  of  time. 

If,  however,  one  is  dealing  with  processes  in  which  wavelengths  less  than 
both  the  Debye  length  and  the  mean  free  path  are  important,  the  path  C cannot  be 
deformed  into  the  path  along  the  real  axis;  this  is  due  to  the  fact  that  the  damp- 
ing frequency  is  higher  than  the  collision  frequency.  The  distribution  function 
at  long  times  will  vary  as  e ^ e *put  while  the  density  will  have  the  more 
rapid  decay  given  by  the  root  tr^  of  the  Eq.  (45).  Thus  the  distribution  function 
for  long  times  contains,  in  addition  to  the  organized  component  of  velocity,  "free- 
particle"  components  which  give  rise  to  negligible  macroscopic  density  fluctu- 
ations. 


For  the  "collision-damping"  case  one  can  easily  write  explicit  expressions 
for  the  asymptotic  behavior  of  the  oscillating  quantities.  We  have 

V = VQ  Re  e ^pX  = VQ  e * cos  (px  - wft)  (47) 


4iren  w.t 

E = - e sin  (px  - <a  t) 

x p r ' 


(48) 


<t>  = 


V e r r uu  uw 

~ 7” ~l  cos  (px-«rt)  i(A.+wi)(l+  £)  X -(pu-«r)(l-— ^i)  — P l 

w.)  +(pu-wr)  I l pa  Up  pa  J 

r . UW  . . kw..  UW^  1 ”| 

» (px-wft)  A(pu-«r)  (1+ 2)  -(^+«i)(1-— zj f \ 

I pa  w pa  J J 


(49) 
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The  real  and  imaginary  parts  of  u are  given  by  the  equations  obtained  in  our 
study  of  the  dispersion  relation  in  Secs.  5 to  7. 

Using  the  above  expressions  it  is  now  easy  to  verify  that  the  approxi- 
mations made  in  the  linearization  procedure  are  valid  for  all  velocities  in  this 
long  wavelength  case.  The  amplitude  VQ  can  be  chosen  arbitrarily  small  and, 
provided  the  collision  frequency  is  not  zero,  the  linearization  leads  to  consistent 


results . 


1 C C 

For  the  high-pressure  limit  we  have  from  Eq.  (68)  t*>  + ^ ^ ; 


u.->»0  as  X.-^oo.  The  distribution  function  is 
i 


♦ Vo  COS  (pX  “ V*  {l  + "2  ^ 


3 2 2 ^ 2 2 

For  the  low-pressure  plasma  u u>  + -«•  P--v  ; w.  -»•  - — P + as  X.  — 0. 
r rp22i  oj  2 


Then 


P <*> 


I 

vo e 1 r 

^(pu-Wj?  I 


cos(px-cj  t)  — P (pu-oa  ) - sin(px 
r l r 

pa 


•V) 2 — # 

pa  J 


For  sound  waves  in  an  uncharged  gas  we  again  have  4*  — ► VQ  cos  (px-«rt)^l+—  ^ 

c.  r a4  P 

here  however,  one  must  take  u>r  = pa.  Thus 

<t>  VQ  cos  p (x  - at)  (1  + £)  . (52) 

All  the  above  distribution  functions  exhibit  the  same  general  feature  that  the 
magnitude  is  greater  for  those  particles  moving  in  the  direction  of  the  wave  than 
it  is  for  those  moving  in  the  opposite  direction.  This  is  in  fact  necessary  for 
maintaining  a wave.  For  the  low-pressure  plasma  case  we  note  also  the  oppo- 
site contributions  from  particles  moving  faster  and  slower  than  the  wave. 


5.  DISPERSION  RELATION  FOR  INITIAL  VALUE  PROBLEM 

The  present  section  is  devoted  to  a discussion  of  the  techniques  needed 
for  the  study  of  the  dispersion  relation  Eq.  (45).  The  results  for  important 
special  cases  are  given  in  the  following  two  sections.  At  this  stage,  it  is  con- 
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venient  to  introduce  dimensionless  variables.  We  measure  frequencies  in  units 
of  the  plasma  frequency  co^,  velocities  in  units  of  the  thermal  velocity  a = (kTQ/m)^^, 
and  lengths  in  units  of  the  Debye  length  = (kT/4TrnQe^)^^.  In  these  units  the 
dispersion  relation  becomes  2/^ 

ip  = i-y-  J + ipk  - <rX.C  1 d£  . (53) 

(2tt)1/Z  Jc  \ + <r  + ip£  L J 

We  shall  at  times  deal  with  w rather  than  <r;  the  two  quantities  are  connected  by 


<r  = - 10). 


Let  us  first  study  the  position  of  the  singularity  of  the  integrand  in  the  b 
plane.  Since  the  oscillations  are  damped,  <r  has  a negative  real  part  and  a neg- 


ative imaginary  part.  Thus  X.  + <r  can  be  either  positive  or  negative  according 
as  the  damping  frequency  is  less  than  or  greater  than  the  collision  frequency. 
Since  p is  real  this  means  that  the  singularity  C.  - U + «r)/ip  is  in  either  the 
first  or  the  fourth  quadrant  in  the  C plane.  For  in  the  first  quadrant,  the 

contour  appropriate  to  the  e ^pX  substitution,  i.e.,  along  the  real  axis. 


gives  correct  results.  This  represents  the  situation  where  the  main  damping 
arises  from  collisions.  When  CQ  is  in  the  fourth  quandrant,  the  damping 
frequency  is  larger  than  the  collision  frequency;  the  additional  damping  is  of  the 
Landau  type  and  results  from  the  destruction  of  phase  by  drift  motion.  The 


representation  of  this  latter  type  of  damping  requires  the  use  of  the  contour  C. 
The  damping  frequency  is  small  for  long  wavelengths  but  increases  rapidly  as 
the  wavelength  approaches  the  Debye  length. 

The  most  direct  approach  to  the  evaluation  of  Eq.  (53)  would  appear  to  be 
an  expansion  of  the  term  l/(X.  + ar  + ipC  ) in  powers  of  ip  C /(X.  + tr)  and  an  inte- 
gration term  by  term.  This  is  a development  in  powers  of  p/(X  + c)  or,  revert- 
ing to  usual  units,  in  powers  of  pa/(X  - iw).  Thus  one  requires  pa  to  be  small 
compared  to  either  the  collision  frequency  or  the  frequency  of  the  waves.  Alter- 
natively, the  wavelength  is  long  compared  to  the  Debye  length  or  the  mean  free 


- 28  - 


path.  This  is  actually  the  case  in  two  important  limits.  The  case  of  sound 
oscillations  in  neutral  gases  is  one  where  the  collision  frequency  is  large 
compared  to  pa  for  all  but  very  short  waves  and  very  low  pressures.  A second 
situation  occurs  with  electron  plasma  oscillations  of  wavelength  long  compared 
to  the  Debye  length  ( ~ £^.  = | ^ ng  * ) irrespective  of  the  col- 

lision time,  so  that  one  can  study  the  entire  range  of  pressures.  In  fact,  both 
of  these  cases  can  be  treated  by  the  same  expansion.  The  expansion  breaks 
down  for  heavily  damped  waves  in  both  the  neutral  and  ionized  gas  cases.  An 
objection  to  the  above  procedure  is  that  the  series  does  not  converge;  it  is,  in 
fact,  asymptotic  so  that  it  is  useful  where  the  first  few  terms  are  all  important. 
In  order  to  establish  this  and  to  give  alternative  developments,  we  express  the 
dispersion  relation  in  terms  of  the  error  function  of  complex  argument. 

Let  us  now  displace  the  contour  parallel  to  the  real  axis  so  that  it  has  an 
imaginary  co-ordinate  equal  to  that  of  the  singularity.  If  one  breaks  up  the 
integral  into  the  contribution  of  the  small  semicircle  about  the  singularity  and 
that  of  the  principal  value  of  the  integral  one  finds 


Making  use  of  the  identity 


r 

V n 


V'2  si 


sinh  9 


£ 


and  introducing 


r\  i0  _ k - i« 


we  find  the  dispersion  relation 

' 1 + i<A  = 1 + 


_-a‘2  .rvl  a2 

Jo 


* 

e d Qt  , e 


(56) 

(57) 

(58) 


Here 


F (&)  = i 


(59) 
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is  the  error  function  of  complex  argument. 

Examination  of  Qi  shows  that,  for  waves  with  positive  real  frequency, 
it  lies  either  in  the  fourth  or  third  quadrant.  The  function  F (^2)  is  treated 
in  detail  by  Rosser.  has  the  expansion 

£l2  f 


F (X2) 


„ _ 2n  + 1 

2 (-i)n  £ 

n = 0 (2n  + l)n  j 


(60) 


convergent  for  all  values  of  £2 , but  useful  only  for  values  of  X2  with  modulus 
less  than  unity.  F (i2)  can  also  be  expanded  in  an  asymptotic  series  of  the  form 


F (£2)  = - —U,  + — , - 1-‘-3  ' y + . . . 

z£l  8 X25  16  ft7 


(61) 


This  asymptotic  series  is  in  powers  of  E a -^i2"  and  yields  the  same  result  as 
the  procedure  outlined  at  the  beginning  of  this  section. 

Finally,  we  note  that  F (X2)  can  also  be  expanded  as  a continued  fraction 

f (fi)  = -W  4r  -?r  • • • (62) 

Z\l  + 2 M + 2ii+  2 

for  & in  the  first  and  fourth  quadrants.  This  expansion  does  not  correspond  to 
a simple  development  in  powers  of  the  Debye  length  or  mean  free  path  divided 
by  the  wavelength;  however,  it  converges  rapidly  over  a large  portion  of  the 
complex  ^2  plane  and  is  more  generally  useful  than  the  asymptotic  development. 
To  obtain  information  in  the  regions  of  moderate  and  strong  absorption  for  both 
neutral  and  ionized  gases,  we  require  a knowledge  of  F (ii>  over  portions  of  the 
complex  £2  plane  not  covered  by  any  of  the  developments  given;  this  knowledge 
has  been  obtained  by  numerical  computations. 

We  also  require  a knowledge  of  F<&)  for  £2  in  the  second  and  third 
quadrants.  This  is  found  by  making  use  of  the  identity 

a2 


F (-£2)  = IT1/2 


f (X2) 


(63) 


16)  J.  B.  Rosser,  "Theory  and  Application  of  J 
Brooklyn,  N.Y.,  1948. 


-x 


dx,  " Mapleton  House, 
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We  see  that  the  damping  frequency  is  proportional  to  p , in  contrast  to  the 

damping  arising  from  the  electron-neutral  gas  collisions,  which  contains  a 

leading  term  independent  of  wavelength.  The  coefficient  of  p a /u‘  is  plotted 

in  Figs.  1 and  2.  The  absorption  reaches  a maximum  when  X.  = Up,  i.e.,  when 

the  collision  frequency  equals  the  plasma  frequency.  At  this  point  the  real 

part  of  the  frequency  is  exactly  intermediate  between  the  limiting  values  for 

zero  and  infinite  collision  times.  The  absorption  is  approximately  linear  with 

2 2, 

\ at  the  low-pressure  limit,  varying  as  — -y—  — , while  at  high  pressures  it 

up 

2 2 u p r 

is  hyperbolic,  varying  as  E * — -r-E-  . The  negative  frequency  root  of  Eq.  (67) 


yields 


u = - I u + 


2 2 
pa 


t?  { 


+ 3u^  - i 2X.u 
P P 


and  gives  absorption  and  dispersion  which  corresponds  completely  to  the  results 

for  the  positive  frequency  case-  It  is  important  to  note  that  the  absolute  magni- 

2 2 

tude  of  the  damping  is  very  small  since  it  is  proportional  to  ^ ^ . For  example, 

“ -4 

for  pa  — Wp^QQ  the  maximum  damping  frequency  is  less  than  F 10 

The  dispersion  relation  Eq.  (68)  gives  adequate  information  for  pa  « u^. 

We  are  also  interested  in  the  case  where  pa  is  approximately  or  greater  than 

Up,  i.e.,  waves  of  length  comparable  to  the  Debye  length.  For  this  case  the 

development  in  powers  of  pa/up  is  not  appropriate;  one  must  also  have  detailed 

numerical  information  concerning  F (i2).  In  Figs.  3 and  4 we  plot  the  real  and 

imaginary  parts  of  u for  the  case  pa  = Up  (waves  equal  to  the  Debye  length)  as  a 

function  of  the  collision  frequency.  The  results  are  the  same  as  for  the  long- 

wavelength  case  when  X.  » u . However,  at  values  of  X.  smaller  than  u the 

P P 

very  large  "drift"  damping  discussed  in  Sec.  4 manifests  itself. 
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7.  OSCILLATIONS  OF  AN  UNCHARGED  GAS 


The  general  results  of  Secs.  3 to  5 apply  to  the  initial-value  problem 
for  uncharged  gases;  the  plasma  frequency  is  then  zero.  For  our  constant- 
collision  time,  isothermal  model  the  general  dispersion  relation  Eq.  (58)  be- 


comes, in  ordinary  units, 

2 2 

- P- 4-  = l + 

1WA. 


X - iw 


pa  (2) 


As  discussed  in  Sec.  5 the  asymptotic  expansion  Eq.  (61)  is  valid  provided 
pa  « X,  i.e.,  the  frequency  of  the  sound  wave  is  smaller  than  the  collision 
frequency.  The  first  few  terms  of  this  expansion  yield  the  result 

i/  + p2a2  f . 1 + 4iuA  1 + P4a4  f 3 - 18i“A  1 = 0 . (70) 

j (7  - WM*  ] (1  - WM4  I 

2 2 2 

In  the  limit  of  infinite  collision  frequency  we  find  u>  = p a , the  correct  isother- 
mal sound  speed.  The  correction  for  finite,  but  large  collision  frequency,  ex- 


pressed in  powers  of  pa/X,  is 


w = pa 


This  expression  indicates  both  an  absorption  and  dispersion  of  sound  for  finite 
value  of  the  collision  frequency.  The  physical  situation  is  similar  to  the  Lang- 
muir description  of  the  damping  of  plasma  waves  with  wavelength  near  to  the 
8) 

Debye  length.  ' Here,  because  of  the  translational  motion,  fast  moving  par- 
ticles can  occasionally  move  ahead  of  the  sound  wave  from  one  region  to  a 
distant  region  without  undergoing  a collision  and  thus  carry  ordered  motion  from 
one  region  to  another  where  the  phase  for  ordered  motion  is  different.  The 
ordered  energy  is  then  dissipated  into  random  thermal  energy  by  collisions. 


This  effect,  of  course,  increases  with  increasing  mean  free  path  and  with  de- 
creasing wavelength. 
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The  series  expansion  Eq.  (70)  breaks  down  when  pa  — X.  and  one  must 
resort  to  numerical  evaluation  of  Eq.  (69).  The  results  for  the  real  and  imagi- 
nary parts  are  presented  in  Fig.  5.  The  accuracy  of  these  results  is,  however, 
limited  by  the  inadequate  tables  of  the  function  F (£1)  available  to  us  at  present. 

According  to  the  discussion  of  Sec.  4 the  asymptotic  behavior  of  the 
distribution  function  differs  according  as  = - ia»o  = - k - ipu,  or  = - iw^  , 
as  given  by  the  dispersion  relation  Eq.  (69),  has  the  larger  real  part.  The 
second  case  is  what  we  have  previously  called  the  "collision-damping"  case. 

For  this  case  the  damping  frequency  is  less  than  the  collision  frequency  and  the 
asymtotic  behavior  of  the  distribution  function  is  the  same  as  that  of  the  density. 
From  Fig.  5 we  see  that  this  is  the  situation  when  pa  < 1.  5X..  For  the  first 
case  the  damping  frequency  is  greater  than  the  collision  frequency,  the  distri- 
bution function  and  density  have  different  asymptotic  behaviors,  and  the  distri- 
bution function  is  directly  sensitive  to  the  initial  value  of  the  distribution  function. 
The  individual  features  of  the  particular  initial  value  or  boundary- value  problem 
treated  are  then  important. 

For  the  "collision-damping"  case,  one  may  study  the  physical  processes 
with  the  aid  of  the  distribution  function  Eq.  (49).  Using  the  value  of  (■>  given  by 
Eq.  (71)  we  see  that  the  deviation  from  the  Maxwellian  distribution  is  greatest 
for  fast  moving  particles  and  short  wavelengths.  From  Eq.  (39)  for  the  mass 
velocity,  we  see  that  the  greatest  contribution  to  the  out-of-phase  component  of 
mass  velocity  comes  from  fast  moving  particles  and  short  wavelengths  in 
accordance  with  the  picture  discussed  above. 

8.  GENERAL,  TREATMENT 

In  this  section  we  shall  remove  the  restriction  made  previously  and  allow 
for  temperature  fluctuations.  The  problem  is  now  characterized  by  the  full  set 
of  Eqs.  (15)  to  (20).  The  passage  to  the  linear  approximation  involves  the 
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additional  fluctuating  quantity  T(x,  t),  defined  by 

T (x,  t)  = Tq  |l  + T ( x,  t)|  (72; 

Equations  (21)  to  (25)  remain  valid  in  the  general  case,  but  the  expansion  of  <f> 


in  Eq.  (26)  must  be  replaced  by 


$>  = F { 


v * q + T 


(#  -1)} 


Equation  (19)  defining  the  temperature  leads  to 

3kT 

—2-(1/  + T)=/vZ  ?4>dv  • (74) 

The  basic  set  of  linear  equations  governing  the  Fourier  transformed  vari- 


ables for  a disturbance  in  the  x direction  is 


e(E  ) 

' p'x 


» x r 


f - 4>  + -T-7 

P P k- 


( V*  ° / “ ♦p ? d 


— ( V + T ) =/v^  <j>  F dv 
P P J P 

-*  + -* 

ip(E  ) = - 4iren  V 

r ' p'x  o p 


One  may  eliminate  (E  )x  and  (q^)x  in  the  same  way  as  in  the  isothermal  case. 

Corresponding  to  Eq.  (44)  one  now  finds  a coupled  set  of  equations  for  the 

Laplace  transforms  of  V and  T . This  set  is 

P P 
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f V2  F dv  r 

. Xm 

— + 

,2  1 
k u 
0 

/ X + <r  + ipu  L 

* " FT 

■ o 

ip 

iP  J 

, _ f , Xm  f F d v f mv^  3 \ ] 

+ Vo-  ) 1 " iFr  Jr+TTipu  lurr " z)  j 

V-  ..  /n\ 

} 


(77) 


m 

TFT 

o 


2 - ♦ r u 
v F dv  j Xm  ^ 


[ v F 

J X + cr  + 


t et;  -tp~  + V 0) 


The  dispersion  relation  in  dimensionless  units  is  then 

{ 


In  order  to  study  this  dispersion  relation  we  express  all  the  integrals 
occuring  in  terms  of  the  single  integral 

r -£2/2  t 

1 J ^ 


mP1 

which  was  discussed  in  Sec.  5. 


i - 9 


(79) 
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^ ///  ^ ifc-V2  - *fl4> 

J 


(2tr) 

We  also  find 


_1  fff  «'Y, /.I -d-?  |\  + *22  + 4-1  = - i^LzJ^flEk 

(2ir)3/2  JJJ  v + i(Pu  - «)  |_  P ipj  p2 

+ F(fl)  j^I+Ap(l  + X2-p>l/?fl)  | , 

-V2/2  2 

• \T 


J 


- (1  + X2  - p\/?Xl)  [3  - 2^i  - 4X2(1  -X2)f] 

P 

Using  these  results  we  can  write  the  dispersion  relation  as 

| ! + F (ft)  [^1  + i^u  + x2  - px/zfl)] 

|l  + £ ^ [ - 2X2  + 2&  +(-2  + 2X1  - 4X2  ) F ]J. 


(80) 


(81) 


{ 


1 - ^ [ XI  + 2(1  - XI  )f]  + ~^~2  (1  + X.2  - p\Jz  XI)  [ 3 - 2Xi  - 4X2(1-  XI 

P 3p 


| I f [ -U  + (1  + zXl2)  F]  j 


The  dispersion  relation  Eq.  (81)  for  the  general  case  is  considerably  more 
involved  than  the  relation  Eq.  (58),  which  holds  in  the  isothermal  case.  However, 
as  emphasized  in  Sec.  5 one  can  obtain  simple  results  for  wavelengths  long  com- 
pared to  either  the  Debye  length  or  the  mean  free  path.  For  this  case  the  asymp- 
totic development  for  F (X2)  holds  and  one  finds  (in  dimensionless  units) 


(!-«)(-  iu)  + 


_M>2  f 


5 19w  5iu>X  , 3iw 


iu>  1 

x J 


( X.  - iw) 


{■ 


5A.3  - iw\2  + X(7  + 3w2)  - 20(1  - w2)  - 15iw  l 


J 


(82) 
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For  the  study  of  the  oscillations  of  an  ionized  gas  we  obtain  an  expression  for 

2 2 2 

the  frequency  which  corresponds  to  Eq.  (66).  Keeping  terms  of  order  p a 
we  find  for  the  frequency 


2 2 , 

, p a 1 

W = « + ■=* —x T 

p Zw  /.  2 2. 

p (X  + u>  ) 


XT  + 3 co 


• 4 ^ 

” 1 **■ 

3 p 


} 


(83) 


At  high  pressures  the  collision  frequency  X is  large  compared  to  the  plasma 

2 2 
n ^ 

frequency  and  one  finds  co  co  + . This  is  the  value  found  previously 

^ P 11) 

by  a generalization  of  the  Thomson  method  ' and  verifies  the  high-pressure 

nature  of  the  transport  type  of  procedure.  At  low  pressures  we  find  co  co^  + 

2 2 
P a 

2 co 


3,  the  correct  dispersion  relation  for  waves  of  length  large  compared  to 


the  Debye  length.  The  maximum  damping  frequency  for  the  general  case  occurs 
again  when  the  collision  frequency  is  equal  to  the  plasma  frequency.  It  is  in 
magnitude  2/3  that  for  the  isothermal  case;  the  decreased  absorption  accompanies 
the  smaller  dispersion  for  the  general  case.  The  detailed  variation  of  frequency 
as  a function  of  collision  frequency  as  given  by  Eq.  (83)  is  plotted  in  Figs.  1 and  2. 
In  the  event  that  the  wavelength  is  comparable  to  both  the  collision  frequency  and 
the  Debye  length  one  must  resort  to  numerical  information  concerning 
Tables  of  sufficient  accuracy  to  solve  the  complicated  Eq.  (81)  are  not  available 
to  us  at  present;  we  leave  the  detailed  investigation  to  a later  time. 

For  an  uncharged  gas  the  expression  Eq.  (82)  reduces  to 


W2  + 


2 2 

-EJL_ 


fl9  iu 

- 

4 4 

r + P a 

f 5iX  20 

3iw  | 

lT  X 

\2(1  - iw/X)4 

T " 

(1  - iu/X) 

Developing  in  powers  of  pa/X,  one  finds  for  pa/X  <<  1 

!/2  f /svl/2 


(84) 


w=(!)  pa - 1 (4) 


P-  I + U.3  J2L*! 

X 5 *2 


} 


(85) 


Thus  one  finds  both  an  absorption  and  a dispersion  of  sound  waves.  No  attempt 
has  been  made  in  this  paper  to  connect  our  collision  frequency  X with  molecular 
parameters;  it  is  known  therefore  only  to  within  a factor  of  the  order  of  unity. 
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This  comparison  will  be  undertaken  in  Paper  4,  where  we  also  discuss  models 

more  realistic  than  the  present  constant-collision-time  model.  For  the  present 

we  note  that  if  one  fits  X.  by  taking  the  absorption  equal  to  that  predicted  by  the 

90  a^ 

Navier-Stokes  equations  (and  also  the  Burnett  equations)  one  finds  X = — , 

1 / 8kT \ ^ ^ * a ' 

where  7J  = -g  PQ  — j / ; 1)  is  the  coefficient  of  viscosity;  Z is  the  mean 

free  path.  The  coefficient  of  the  dispersive  term  in  Eq.  (85)  is  then  25  percent 
higher  than  that  given  by  the  Navier-Stokes  equation;  it  thus  lies  between  the 
predictions  of  the  Navier-Stokes  and  of  the  Burnett  equations.  The  interesting 
experimental  observations  of  Greenspan  on  the  dispersion  and  absorption  of 
helium  gas  reveal  important  discrepancies  between  experiment  and  the  predictions 
of  both  the  Navier-Stokes  and  the  Burnett  equations  when  pa  — X.  From  the  point 
of  view  of  the  present  work  the  study  of  this  region  requires  a numerical  analysis 
of  Eq.  (81).  We  shall  therefore  later  take  up  the  study  of  the  translational  dis- 
persion in  a more  detailed  way. 


